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Abstract: We present a numerical study of the ghost-gluon vertex and of the correspond- 
ing renormahzation function Zi(p^) in minimal Landau gauge for SU{2) lattice gauge the- 
ory. Data were obtained for three different lattice volumes (V = 4'^,8^,16^) and for three 
lattice couplings (3 = 2.2,2.3,2.4. Gribov-copy effects have been analyzed using the so- 
called smeared gauge fixing. We also consider two different sets of momenta (orbits) in 
order to check for possible effects due to the breaking of rotational symmetry. The vertex 
has been evaluated at the asymmetric point (0; p, —p) in momentum-subtraction scheme. 
We find that ^i(p^) is approximately constant and equal to 1, at least for momenta p > 1 
GeV. This constitutes a nonperturbative verification of the so-called nonrenormalization of 
the Landau ghost-gluon vertex. Finally, we use our data to evaluate the running coupling 
constant as{p^). 
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1. Introduction 

In the framework of quantum field theory, Faddeev-Popov ghosts are introduced in order to 
quantize non-Abelian gauge theories. Although the ghosts are a mathematical artifact and 
are absent from the physical spectrum, one can use the ghost-gluon vertex and the ghost 
propagator to calculate physical observables, such as the QCD running coupling as(p^), 
using the relation [1] 

„,(,2)=„„M£)a«. (1.1, 

Here a.o = gfJ^iT is the bare coupling constant and Z^{p'^), Z'^{p'^) and Zi{p^) are, respec- 
tively, the gluon, ghost and ghost-gluon vertex renormalization functions. In Landau gauge, 
the vertex renormalization function Ziijp') is finite and constant, i.e. independent of the 
renormalization scale p, to all orders of perturbation theory (see [2] and the nonrenormal- 
ization theorems in [3, Chapter 6]). However, a direct nonperturbative verification of this 
result is still lacking. This is the purpose of the present work. Let us stress that, if Ziip^) is 
finite and constant also at the nonperturbative level, the equation above can be simplified, 
yielding a nonperturbative definition of the running coupling constant that requires only 
the calculation of gluon and ghost propagators [4, 5]. A nonperturbative investigation of the 
structure of the ghost-gluon vertex is also important for studies of gluon and ghost propa- 
gators using Dyson-Schwinger equations (DSE) [6]. In fact, in these studies, one makes use 
of Ansatze for the behavior of the propagators and vertices in the equations, in order to 
obtain solvable truncation schemes. Of course, a nonperturbative input for these quantities 
is important for a truly nonperturbative solution of the DSE. 
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Let us recall that in Landau gauge the gluon and ghost propagators can be expressed 
(in momentum space) as 

= S''' {S,. - ^) D{q') = 5''^ [5,. - ^) ^ (1-2) 

G''%q-q) = -<5^^ G{q^) = -5^^ ^ , (L3) 

where F{q'^) and J{q^) are, respectively, the gluon and ghost form factor and the color 
indices b and c take values 1,2, .. . , A''^ — 1 in the SU{Nc) case. Then, in tlie momentum- 
subtraction scheme one has that the gluon and ghost renormalization functions are given 

by 

^ z^\p^)F{q^) (L4) 

jRiq\p'') = Z^\p'')J{q'') (L5) 

with the renormalization conditions 

FRip^p^) = Jr{p\p^) = 1 . (1.6) 

Thus, if one can set Zi{p'^) = 1, the running coupling (1.1) can be written as [4, 5] 

as{p^) = aoF{p^)j\p^) , (1.7) 

i.e. one only needs to evaluate the gluon and ghost form factors defined above. In recent 
years, the infrared behavior of these form factors has been extensively studied (in Landau 
gauge) using different analytical approaches [4, 5, 7, 8, 10, 9, 11, 12, 13, 14, 15, 16, 17, 18, 
19, 20, 21, 22, 23, 24, 25]. Also, numerical studies of these form factors and of the running 
coupling defined in eq. (1.7) have been reported in [26, 27, 28, 29] for the SU{2) group and 
in [30, 31, 32, 33, 35, 36, 34, 37] for the ^[/(S) case. 

An indirect evaluation of the ghost-gluon vertex renormalization function has been 
recently presented in [29], confirming that ZiijP') is finite in the continuum limit. On the 
other hand, a direct nonperturbative verification of this result from a numerical evaluation 
of the ghost-gluon vertex is still missing. Let us stress that a direct evaluation of Ziip^) 
would allow a study of the running coupling constant [38] using eq. (1.1) instead of eq. 
(1.7). Such a study may improve the precision of the determination of agip^), since in that 
case one does not need, in principle, to use the so-called matching rescaling technique [29] 
when considering data obtained at different /? values. 

In this paper we study numerically the ghost-gluon vertex and the corresponding renor- 
malization function Zi(j)'^) for the SU{2) case in the minimal Landau gauge. The definition 
of this renormalization function (in the continuum) is presented in Section 2. Numerical 
simulations are explained in Section 3. In particular, in Section 3.1 we define the vertex 
renormalization function on the lattice and compare our direct determination of Zi(p^) to 
the indirect evaluation presented in [29]. Note that, since the numerical evaluation of the 
ghost-gluon vertex may be affected by Gribov-copy effects (see discussion in Sections 3 and 
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Figure 1: Notation (momenta and indices) for the ghost-gluon vertex. 



4 below), we evalute this quantity using two different gauge- fixing methods, leading to two 
different sets of Gribov copies. A comparison of the results obtained in the two cases allows 
us to estimate the influence of Gribov copies (the so-called Gribov noise) on the considered 
quantity. 

Results for the ghost-gluon vertex renormalization function, the gluon and ghost prop- 
agators and the running coupling constant Q!s(p^) are reported in Section 4. Finally, in 
Section 5 we draw our conclusions. Preliminary results have been reported in [39]. 

2. The ghost-gluon vertex 

Following the notation in Ref. [6] (see also Fig. 1) the 3-point function for A'^ (gluon), rj^ 
(ghost) and f]'^ (anti-ghost) fields — corresponding to the ghost-gluon vertex — is given by 

F;^^(x,y,z) = {A'^^{x)ri\y)fj''{z)) . (2.1) 

Going to momentum space and using translational invariance for the 3-point function, one 
gets 

= j d^Z 6-^^=+?-^)^ j d^xdS e.-'^kx+qy) ^ ^a^^^ ^b^y^ -c(q) ^ (2.3) 

= i27r)U\k + q-s)G^'^'^ik,q) , (2.4) 

where S^{k + q — s) implies conservation of momentum. Then, the ghost-gluon vertex 
function is obtained by "amputating" the corresponding 3-point function (see Fig. 2) 



G^'%k,q) 



where s = k + q and the functions D{k'^) and G{q^) have been defined in eqs. (1.2) and 
(1.3), respectively. At tree level (in the continuum) one obtains [1] 

Tt'^iKq) = igor'"'^. (2.6) 
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Figure 2: The 3-point function G^^{k, q) and its relation with the full ghost-gluon vertex T'^'^ik, q). 



where ^'o is the bare couphng and /"^'^ are the structure functions of the SU{Nc) Lie 
algebra. This implies the well-known result that the ghost-gluon vertex is proportional to 
the momentum q of the outgoing ghost. More generally, we can write the relation [6] 



abc 1 



(2.7) 



where Tfj_{k,q) is the so-called reduced vertex function. Then, multiplying the previous 
equation by f<^^<=^S'^^ summing over a,b,c,d and using the relation ^^^^^ jdbcjabc _ ]\f^§da 
we get 



a,b,c 



Since at tree level one has r^(fc, = we can also write 



(2.8) 



(2.9) 



and it is the scalar function r{k^,q^) that gets renormalized, when considering a given 
renormalization scheme. 

Clearly, from the above relations we obtain 



T{k',q') = — ^g^r^(A;,g) = 



—I 



goN,{N^ - 1) q 



(2-10) 



a,6,c 



In our simulations we considered the asymmetric point with zero momentum for the 
gluon (A; = 0), implying s = q. Formula (2.10) then becomes 



■iq') = r(o,g2 



5oiVc(iV2 - 1) q 



(2.11) 



a,b,c 



Let us note that the factor —i in the equations above disappears, since only the imaginary 
part of the vertex function r^''^(0, q) contributes to the quantity '^{q'^) , i.e. we can write 
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Figure 3: Relation between the full (left) and the bare (right) vertices. 



Finally, in momentum-subtraction scheme one fixes the vertex renormalization function 
^lip^) by requiring (see also Fig. 3) 

Eniq^p') = Z,{p')^i<l') (2.13) 
= 1 , (2.14) 

namely the renormalized reduced vertex function r^(0, g) is equal to the tree-level value 
ig^ at the renormalization scale. Then, the vertex renormalization function is given by 

Z^\p') = E{p') . (2.15) 

We also obtain that the running coupling can be written as 

gsip^) = go Z^\p'') Zy\p^) Up'') (2.16) 
= 5oH(p2)Fi/2(p2)j(p2) ^2.17) 

where we used eq. (1.3). 

3. Numerical Simulations 

For the simulations (see [40] for details) we consider the standard Wilson action, thermal- 
ized using heat-bath [41] accelerated by hybrid overrelaxation [42, 43, 44]. Since we are 
considering the SU{2) gauge group we parametrize the link variables Ufj,{x) as 

U„{x) = ul{x)l. + iUf,{x)d , (3.1) 

where IL is the 2x2 identity matrix, are the Pauli matrices and the relation 

[ul{x)f + [u,{x)f = 1 (3.2) 
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Lattice Volume V = N'^ 


44 


8^ 


16^ 


No. of Configurations 


1000 


400 


250 



Table 1: Lattice volumes and total number of configurations considered. 



is satisfied. We also consider the standard definition of the lattice gluon field, i.e. 

A,{x) = = EM_EM, (3.3) 

implying that the lattice gluon field component A^^(x) is equal to uj^(x). The connection 
between the lattice link variables U^{x) and the continuum gauge field A^{x) is given by 

U^{x) = e^p[iagoA''^{x)t^'\ , (3.4) 

where = are the generators of the SU(2) algebra, a is the lattice spacing and go 
is the bare couphng constant. Then, the lattice quantity 2^|^(x)/(ag'o) approaches the 
continuum-gluon-field component A^^^ (x) in the naive continuum limit a . Let us also 
recall that a generic lattice momentum q has components (in lattice units) 

q, = 2 sin (^) . (3.5) 

Here = a Nf^^ is the physical size of the lattice in the fi direction, the quantity if^j takes 
values [ — ^ J + 1 ) • • • ) L ^ J ^^cl Nfj^ is the number of lattice points in the n direction. If 
we keep the physical size constant and indicate with q'^ the momentum components in 
the continuum, we find that the lattice components = aq^ take values in the interval 
(— TT, tt] when a — > 0. 

In Table 1 we show the lattice volumes used in the simulations and the corresponding 
numbers of configurations. For each volume we have considered three values of the lattice 
coupling, namely /3 = 2.2, 2.3, 2.4. The corresponding string tensions in lattice units [29] 
are, respectively, equal to a a? = 0.220(9), 0.136(2) and 0.071(1). 

The minimal (lattice) Landau gauge is implemented using the stochastic overrelaxation 
algorithm [45, 46, 47]. We stop the gauge fixing when the average value of (V • A)"^ — see 
Eq. (6.1) in [47] for a definition — is smaller than 10^^^. The gluon and ghost propagators 
in momentum space are evaluated using the relations 



b,fi 

= IJ2 ^ti^ (3-7) 

G{q) = \Y. ' (3.8) 

^ h 

where 

= J]exp[i^(x-j/)](4(x)A^(y)), (3.9) 
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V = 44 


V = 8^ 


V = 16^ 


f3 = 2.2 
l3 = 2.3 
/?= 2.4 


4.3 % 
4.4% 
4.1 % 


25.0 % 
16.8 % 
9.8 % 


97.2 % 
72.8 % 
50.8 % 



Table 2: Percentage of configurations for wiiidi tiie smeared gauge fixing has found a (different) 
Gribov copy. Results are reported for the three lattice volumes and the three /3 values considered. 

G^^(5) = {{M-')'^iq)) , (3.10) 

{M-'f^iq) = ^Yl exp[-ig(a;-y)] {W)^ , (3.11) 
x,y 

V is the lattice volume and M^y is a lattice discretization of the Faddeev-Popov operator 
[{—d + A) ■ d] . The expression of this matrix in terms of the gauged-fixed link variables 
can be found in [48, eq. (B.18)] for the generic SU{Nc) case and in [49, eq. (11)] for the 
SU{2) case considered in the present work. Let us recall that this matrix has a trivial null 
eigenvalue corresponding to a constant eigenvector. Thus, one can evaluate the inverse of 
M^y only in the space orthogonal to constant vectors, i.e. at non-zero momentum. For the 
inversion we used a conjugate gradient method with even/odd preconditioning. 

Finally, in order to check for possible Gribov-copy effects we have done (for each ther- 
malized configuration) a second gauge fixing using the smearing method introduced in [50]. 
To this end we have applied the APE smearing process 

U^{x) ^ (1 - w)U^ix) + w^lix) , (3.12) 

where S^(a;) represents the sum over the connecting staples for the link Un{x), followed by 

a reunitarization of the link matrix. We have chosen w = 10/6 and stopped the smearing 

when the condition „ 

Tr 

^ VFi.i > 0.995 (3.13) 

was satisfied. (This is usually achieved with a few APE-smearing sweeps.) Here Wi^i is the 
1x1 loop and the average is taken over all Wi^i loops of a given configuration. The gauge 
fixing for the smeared configuration as well as the final gauge-fixing step have been done 
again using the stochastic overrelaxation algorithm. The smeared gauge-fixing method is 
supposed to find a unique Gribov copy, even though this copy does not always correspond 
to the absolute minimum of the minimizing functional [50]. In our simulations we found 





V = 4^ 


V = 8^ 


V = 16^ 


/?= 2.2 
/3 = 2.3 
P = 2.4 


46.5 % 

38.6 % 
39.0 % 


52.0 % 
52.2 % 
48.7 % 


63.8 % 
65.4 % 
65.4 % 



Table 3: Percentage of (different) Gribov copies for which the smeared gauge fixing has found a 
smaller value of the minimizing function. Results are reported for the three lattice volumes and the 
three /3 values considered. 
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that the number of different Gribov copies obtained with the smeared gauge fixing increases 
with larger lattice volumes and with smaller /3 values (see Table 2). This is in agreement 
with previous studies [51, 52, 53] and it should be related to an increasing number of local 
minima for the minimizing function when the system is highly disordered. We also found 
that the minimum obtained with the smeared gauge fixing is not always smaller than the 
one obtained without smearing, but at larger lattice volumes the smearing approach seems 
to be more effective in getting closer to the absolute minimum of the minimizing function 
(see Table 3). 

3.1 Ghost-Gluon Vertex on the Lattice 

On the lattice, the definition of the vertex renormalization function is obtained as was 
done in the continuum (see Section 2). The only difference is that, from the weak-coupling 
expansion (or "perturbative" expansion) on the lattice, one obtains that the ghost-gluon 
vertex is given at tree level by [54] 

rf (fc,g) = i5or%cos(^^) , (3.14) 

where s = k + q. Clearly, by taking the formal continuum limit a ^ of the quantity 
V'^'^{k,q)/a one recovers the continuum tree-level result (2.6). Then, in the case A; = 0, 
s = g, eq. (2.15) is still valid on the lattice if one sets 



q 



COS 



\ ) 



with 



= 7^77^ (3.17) 

and 

where is defined in terms of in eq. (3.5). Thus, as in the continuum, we can write 

=® - - 1) ? ^J^, E - C(0, f) . (3.19, 




In order to check for possible effects due to the breaking of rotational symmetry, in our 
simulations we use two types of lattice momenta, i.e. qi = q2 = Q3 = 0, cfi = q and qi = 
Q2 =Q3 = Qa = Q- In the following we will indicate the first type of momenta as asymmetric 
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and the second one as symmetric. Clearly, using the relations Li = L2 = L3 = L4 = L, we 
have _ 

g = 94 = 2 sin (3.20) 



in the asymmetric case and 



q^ = 2sm{^^] (3.21) 



g = 4 sin { ^ ) (3.22) 



in the symmetric one. Then, the quantity E!(5) above becomes 



cos 



in the asymmetric case and 



M a,h,c 



in the symmetric one. Using the trigonometric relation sin(2^) = 2 sin(^) cos(^) these 
formulae can be rewritten as 

m = -^^m. (3.25) 

sin (-T- ) 

where 

^ o,6,c 

in the asymmetric case and 

^® = mN.(Nl-l) E - I i: '■.""(0. 5) (3.27) 

in the symmetric one. 

Finally, using eqs. (2.3), (2.4) and (3.11) we obtain 

Gf\^.q) = V{Al{Q) {M-'t{q)) , (3.28) 

where 

X 

and V is the lattice volume. Then, we can write 

, V{AI{0) (M-i)^'(g)) 



9 



With our notation, in the naive continuum hmit a 0, the lattice quantity Aa!^ D{q)/gQ 
[respectively G{q)] approaches the (continuum) propagator D{q^) [respectively G{q^)]. 
At the same time, we have that 2A1{Q)/{ago) Al{Q) and a? {M-^Y" (q) {M'^f" (q). 
This implies that the (normalized) lattice ghost-gluon vertex gQT^^'^{0,q)/{2a) approaches 
the quantity r^'"^(0, g) in the continuum limit. Also, since q ^ aq when a — > 0, it is easy 
to verify that Z^^^q) ='E{q) is dimensionless. 

As said in the Introduction, in Ref. [29] the ghost-gluon renormalization function ^i(p^) 
was indirectly evaluated. To this end the authors considered eq. (1.1), which can be written 
as 

Z3{p',P)zUp\p) = ^ zHp^p) , (3.31) 

where the (3 dependence is now indicated explicitly. The renormalization functions Z^{p'^,P) 
and Z^{p'^,P) can be evaluated using the matching technique described in detail in [55, Sec. 
III]. Then, one can show [29] that the left-hand side of the above equation rises linearly 
with — In (cra^), where a is the string tension. At the same time, for large enough P values 
one finds that [29] 

— \— oc — \n(aa^) + constant, (3.32) 
ao{P) 

where the constant is /3-independent. It follows that Zf{p'^,P) must be finite in the contin- 
uum limit, if the renormalized coupling asip^) is assumed finite. Of course, this approach 
does not allow either an effective determination of the vertex renormalization function or a 
study of its p dependence. These are the main goals of the present work. 



4. Results 

We have evaluated the reduced ghost-gluon vertex function r^'"^{0,p) at the asymmetric 
point {0;p,—p) using eq. (3.30), in the SU{2) case. (Let us recall that when Nc = 2 the 
structure functions /°*^ are given by the completely anti-symmetric tensor e"^'^.) Results for 
the quantity S(p)i defined in eqs. (3.26) and (3.27), and for Z^^(p) = H(p) [see eq. (3.25)] as 
a function oip = p/a (in physical units) are reported in Figs. 4 and 5. We consider for these 
figures the lattice volume V = 16^, for both asymmetric and symmetric momenta. Error 
bars were obtained using the bootstrap method with 250 samples. (We checked that the 
results do not change when considering 500 samples.) We find that the function Xl(p) has the 
same momentum dependence as the tree-level vertex [i.e. p cos(7rpa/L) ~ sin(27rpa/L)] 
and that (consequently) Zi{p) is approximately constant and equal to 1. 

As explained in Section 3 above, in order to check for Gribov-copy effects we have 
considered two different gauge fixing methods for each thermalized configuration, the second 
of which employs the so-called smeared gauge fixing. Let us recall that Gribov-copy effects 
have been found for the ghost propagator in Landau gauge, at least in the small-momentum 
limit [36, 52, 53, 56]. Recently, such effects were also found for the gluon propagator in the 
infrared region [57]. In Fig. 5 are also shown the data for Z^^{p) obtained for the lattice 
volume V = 16* using the smeared gauge- fixing. The Z^^(j)) data for all lattice volumes 
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3 4 
p [in GeV] 

Figure 4: Results for for the lattice volume V = 16^ as a function of p = p/a in GeV, 

considering asymmetric and symmetric momenta. Error bars were obtained using the bootstrap 
method with 250 samples. The dashed Hne represents the tree-level momentum dependence (~ p) 
of the vertex function in the continuum. 



and (3 values considered in this work are reported in Tables 4 (asymmetric momenta) and 
5 (symmetric momenta) considering the two different gauge-fixing methods (without and 
with smearing). We can clearly see that Z^^{p) decreases as the lattice volume increases 
for a fixed coupling P and that for V = 16^ the dependence of Z^^{'p) on the coupling /3, 
on the type of momenta and on Gribov-copy effects is relatively small. This is also evident 
if we try to fit the data for Z^^{p) to a constant (see Table 6). Nevertheless, it is clear that 
data obtained using the smeared gauge-fixing are generally smaller than those obtained 
without smearing. 

It is difficult to explain the small variations in the results for Z^^{p) when considering 
asymmetric and symmetric momenta or the smeared gauge fixing compared to the standard 
gauge fixing. Indeed, several quantities enter the definition of r^^'^(0,p) [see eq. (3.30] and 
therefore the evaluation of Z^^{p). A more detailed study (see Tables 7, 8 and 9) of the 
gluon propagator at zero momentum D{0), the ghost propagator G{p) for pi = P2 = Ps = 
0,P4 = N/2 and the quantity /"^'^ Im ( A^^iO) (M-^)'"' (p) ) as a function of p = p/a 

does not clarify the situation. In particular, the Gribov-copy effects, if present, are always 
small and within error bars. This is not surprising, since we are considering relatively small 
lattice volumes in the scaling region. We can conclude that the small systematic variations 
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0.8 



• P = 2.2 




▲ p = 2.3 




■ P = 2.4 




O P = 2.2 (Sme) 




A ^ = 1.3 (Sme) 




□ P = 2.4 (Sme) 





p [in GeV] 



0.8 



• P = 2.2 




▲ p = 2.3 




■ P = 2.4 




O P = 2.2 (Sme) 




A ^ = 2.3 (Sme) 




□ ^ = 2A(Sme) 






P [in GeV] 



Figure 5: Results for Z^^{p) for the lattice volume V = 16^ as a function of p = p/a, in GeV, 
considering asymmetric (left) and symmetric (right) momenta. In both cases we show data obtained 
considering the two difTerent gauge-fixing methods (without and with smearing). Error bars were 
obtained using the bootstrap method with 250 samples. 



Q. 
B 
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P = 2.2 


A 


P = 2.3 


□ 


P = 2.4 



o 


p-2.2 


A 


P = 2.3 


□ 


P = 2.4 



p [in GeV] 



3 4 

p [in GeV] 



Figure 6: Results for the running coupHng as{p) as a function p = p/a in GeV, considering 
asymmetric (left) and symmetric (right) momenta. Here we use the standard gauge-fixing method 
(without smearing). Error bars were obtained using the bootstrap method with 250 samples. 



in the results of Zi(p) are probably related to correlations among data at different momenta, 
since they have been obtained using the same set of configurations. 

Finally, we can use our data to evaluate the running coupling constant agip^) using eq. 
(1.1). For Zi{p^) we take the fit to a constant reported in Table 6. The results are shown in 
Fig. 6. We can see that data obtained at different (3 values do not lie on a single curve. This 
is probably related to the fact that our data are not at infinite volume. Indeed, this is a 
well-known effect for the gluon field when evaluated at zero momentum [59, 58, 60, 61, 62]. 
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5. Conclusions 

We have presented the first numerical study of the reduced ghost-gluon vertex function 
r^'"^(0,p) and of the renormalization function Zi{p'^) in minimal Landau gauge. We have 
considered tlie SU{2) case and the asymmetric point (0;p, — p). We found that the vertex 
function has the same momentum dependence of the tree-level vertex and that ^i(p^) is 
approximately constant and equal to 1, at least for momenta p ^ 1 GeV. This is a direct 
nonperturbative verification of the well-known perturbative result that Ziip^) is finite and 
constant to all orders of perturbation theory [2, 3]. In particular, using the result obtained 
at the largest value of (3 considered here (i.e. jd = 2.4) we can write Z^^ip^) = 1.02ly (see 
Table 6), where errors include Gribov-copy effects and discretization errors related to the 
breaking of rotational invariance. 

We are now extending this study, considering larger lattice volumes (and therefore 
smaller momenta), the symmetric point k'^ = = , the 3d case and the SU{Z) gauge 
group. A significant improvement of our results for as{p'^) may be expected by considering 
the symmetric point. 

Acknowledgment s 

This work was supported by Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP) through grants: 00/05047-5 (AC, TM) and 03/00928-1 (AM). Partial support 
from Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq) is also ac- 
knowledged (AC, TM). 

References 

[1] Foundations of Quantum Chromodynamics, T. Muta (World Scientific, Singapore, 1987). 

[2] J. C. Taylor, Nucl. Phys. B 33 (1971) 436. 

[3] 0. Piguet and S. P. Sorella, Lect. Notes Phys. M28 (1995) 1. 

[4] L. von Smekal, R. Alkofer and A. Hauck, Phys. Rev. Lett. 79 (1997) 3591. 

[5] L. von Smekal, A. Hauck and R. Alkofer, Annals Phys. 267 (1998) 1 [Erratum-ibid. 269 
(1998) 182]. 

[6] R. Alkofer and L. von Smekal, Phys. Rept. 353 (2001) 281. 

[7] D. Atkinson and J. C. R. Bloch, Phys. Rev. D 58 (1998) 094036. 

[8] D. Atkinson and J. C. R. Bloch, Mod. Phys. Lett. A 13 (1998) 1055. 

[9] J. C. R. Bloch, Phys. Rev. D 64 (2001) 116011. 

[10] D. Zwanziger, Phys. Rev. D 65 (2002) 094039. 

[11] C. Lerche and L. von Smekal, Phys. Rev. D 65 (2002) 125006. 

[12] C. S. Fischer and R. Alkofer, Phys. Lett. B 536 (2002) 177. 

[13] C. S. Fischer, R. Alkofer and H. Reinhardt, Phys. Rev. D 65 (2002) 094008. 



13 



[14] R. Alkofer, W. Detmold, C. S. Fischer and P. Maris, hep-ph/0309077. 
[15] R. Alkofer, W. Detmold, C. S. Fischer and P. Maris, hep-ph/0309078. 
[16] D. Zwanziger, Phys. Rev. D 67 (2003) 105001. 
[17] J. C. R. Bloch, Few Body Syst. 33 (2003) 111. 

[18] J. M. Pawlowski, D. F. Litim, S. Nedelko and L. von Smekal, hep-th/0312324. 
[19] K. I. Kondo, hep-th/0303251. 

[20] K. I. Kondo, Nucl. Phys. Proc. Suppl. 129 (2004) 715. 

[21] R. F. Sobreiro, S. P. Sorella, D. Dudal and H. Verschelde, Phys. Lett. B 590 (2004) 265. 

[22] R. F. Sobreiro, S. P. Sorella, D. Dudal and H. Verschelde, hep-th/0406161. 

[23] A.C. Aguilar, Ph.D. thesis, IFT-Universidade Estadual Paulista, March 2004 (in Portuguese). 

[24] A. C. Aguilar and A. A. Natale, hep-ph/0405024. 

[25] C. S. Fischer and H. Gies, hep-ph/0408089. 

[26] A. Cucchieri, T. Mendes and D. Zwanziger, Nucl. Phys. Proc. Suppl. 106 (2002) 697. 

[27] J. C. R. Bloch, A. Cucchieri, K. Langfeld and T. Mendes, Nucl. Phys. Proc. Suppl. 119 
(2003) 736. 

[28] K. Langfeld et al., hep-th/0209173, in the Proceedings of Gargnano 2002, Quark 

confinement and the hadron spectrum, edited by N. Brambilla and G. M. Prosperi (World 
Scientific, River Edge, 2003) 297-299. 

[29] J. C. R. Bloch, A. Cucchieri, K. Langfeld and T. Mendes, Nucl. Phys. B 687 (2004) 76. 

[30] D. B. Leinweber, J. L Skuherud, A. G. Wihiams and C. Parrineho [UKQCD Cohaboration], 
Phys. Rev. D 60 (1999) 094507 [Erratum-ibid. D 61 (2000) 079901]. 

[31] D. R. Bonnet et al., Phys. Rev. D 64 (2001) 034501. 

[32] H. Nakajima and S. Furui, Nucl. Phys. Proc. Suppl. 119 (2003) 730. 

[33] H. Nakajima and S. Furui, hep-lat/0303024, in the Proceedings of Nagoya 2002, Strong 
coupling gauge theories and effective field theories, edited by M. Harada et al. (World 
Scientific, River Edge, 2003) 67-73. 

[34] S. Furui and H. Nakajima, hep-lat/0309166. 

[35] S. Furui and H. Nakajima, Phys. Rev. D 69 (2004) 074505. 

[36] H. Nakajima and S. Furui, Nucl. Phys. Proc. Suppl. 129 (2004) 730. 

[37] S. Furui and H. Nakajima, hep-lat/0403021. 

[38] S. Furui and H. Nakajima, hep-lat/0012017, in the Proceedings of Wien 2000, Quark 
confinement and the hadron spectrum, edited by W. Lucha and K. M. Maung (World 
Scientific, Singapore, 2001) 275-278. 

[39] A. Mihara, A. Cucchieri and T. Mendes, hep-lat/0408021, presented at the IX Hadron 
Physics and VII Relativistic Aspects of Nuclear Physics Workshops, Angra dos Reis, Rio de 
Janeiro, Brazil (March 28-April 3, 2004). 



- 14 - 



[40] A. Cucchieri, T. Mendes, G. Travieso and A. R. Taurines, hep-lat/0308005, in the 
Proceedings of the 15th Symposium on Computer Architecture and High Performance 
Computing, Sao Paulo SP, Brazil, 2003, edited by L.M. Sato et al. (IEEE Computer Society 
Press, Los Alamitos CA, 2003) 123-131. 

[41] M. Creutz, Phys. Rev. D 21 (1980) 2308. 

[42] S. L. Adler, Nucl. Phys. Proc. Suppl. 9 (1989) 437. 

[43] F. R. Brown and T. J. Woch, Phys. Rev. Lett. 58, 2394 (1987). 

[44] U. WolfT, Phys. Lett. B 288, 166 (1992). 

[45] A. Cucchieri and T. Mendes, Nucl. Phys. B 471 (1996) 263. 

[46] A. Cucchieri and T. Mendes, Nucl. Phys. Proc. Suppl. 53 (1997) 811. 

[47] A. Cucchieri and T. Mendes, Comput. Phys. Commun. 154 (2003) 1. 

[48] D. Zwanzigcr, Nucl. Phys. B 412 (1994) 657. 

[49] A. Cucchieri, Nucl. Phys. B 521 (1998) 365. 

[50] J. E. Hetrick and P. de Forcrand, Nucl. Phys. Proc. Suppl. 63 (1998) 838. 
[51] P. Marenzoni and P. Rossi, Phys. Lett. B 311 (1993) 219. 
[52] A. Cucchieri, Nucl. Phys. B 508 (1997) 353. 

[53] A. Cucchieri and T. Mendes, Nucl. Phys. Proc. Suppl. 63 (1998) 841. 

[54] H. Kawai, R. Nakayama and K. Seo, Nucl. Phys. B 189 (1981) 40. 

[55] A. Cucchieri, T. Mendes and A. R. Taurines, Phys. Rev. D 67 (2003) 091502. 

[56] T. D. Bakeev, E. M. Ilgenfritz, V. K. Mitrjushkin and M. Mueller-Preussker, Phys. Rev. D 
69 (2004) 074507. 

[57] P. J. Silva and 0. Ohveira, Nucl. Phys. B 690 (2004) 177. 

[58] A. Cucchieri, Phys. Rev. D 60 (1999) 034508. 

[59] A. Cucchieri, Phys. Lett. B 422 (1998) 233. 

[60] A. Cucchieri and D. Zwanziger, Phys. Rev. D 65 (2002) 014001. 

[61] A. Cucchieri and D. Zwanziger, Phys. Lett. B 524 (2002) 123. 

[62] A. Cucchieri, F. Karsch and P. Petreczky, Phys. Rev. D 64 (2001) 036001. 



- 15 - 



p (GeV) 


0.366 


0.718 


1.04 


1.33 


1.56 


1.73 


1.84 


= 4 
N = 4, sme 


- 




- 


1.40(5) 
1.39(4) 


- 




- 


N = 8 
N = 8, sme 




1.08(5) 
1.10(5) 




1.11(4) 
1.12(4) 




1.13(5) 
1.14(4) 




N = 16 
N = 16, sme 


1.05(5) 
0.98(5) 


1.13(5) 
1.06(5) 


1.16(5) 
1.07(5) 


1.14(5) 
1.06(5) 


1.14(5) 
1.05(5) 


1.15(5) 
1.05(5) 


1.15(5) 
1.05(5) 


p (GeV) 


0.466 


0.913 


1.33 


1.69 


1.98 


2.20 


2.34 


= 4 
N = 4, sme 








1.31(4) 
1.30(4) 








N = 8 
N = 8, sme 




1.30(5) 
1.24(6) 




1.22(5) 
1.18(5) 




1.21(5) 
1.18(5) 




= 16 
= 16, sme 


1.01(4) 
1.06(5) 


1.05(4) 
1.09(5) 


1.05(5) 
1.08(5) 


1.04(5) 
1.07(5) 


1.04(5) 
1.07(5) 


1.04(5) 
1.07(5) 


1.03(4) 
1.06(5) 


p (GeV) 


0.644 


1.26 


1.83 


2.34 


2.75 


3.05 


3.24 


A^ = 4 
N = 4, sme 








1.36(4) 
1.35(3) 








AT = 8 
N = 8, sme 




1.16(5) 
1.15(5) 




1.06(5) 
1.07(5) 




1.06(4) 
1.06(4) 




L = 16 
A^ = 16, sme 


1.11(6) 
1.08(6) 


1.10(5) 

1.06(5) 


1.08(5) 
1.04(5) 


1.08(5) 
1.05(5) 


1.07(5) 
1.05(5) 


1.08(5) 
1.05(5) 


1.08(5) 
1.05(5) 



Table 4: Results for Z^'^ip) as a function of p = p/a (in GeV) for the three lattice sides V = 4, 8, 16 
and the three /? values considered (i.e. 2.2, 2.3 and 2.4, respectively top, center and bottom part of 
the table) in the case of asymmetric momenta, considering the two different gauge-fixing methods 
(without and with smearing). Error bars were obtained using the bootstrap method with 250 
samples. 



16 



p (GeV) 


0.732 


1.44 


2.08 


2.65 


3.12 


3.47 


3.68 


N = A 


- 


- 


- 


1.22(4) 


- 


- 


- 


N = A, sme 


- 


- 


- 


1.23(4) 


- 


- 


- 


N = 8 


- 


1.07(5) 


- 


1.01(5) 


- 


0.99(5) 


- 


N = 8, sme 


- 


1.03(5) 


- 


0.98(5) 


- 


0.97(5) 


- 


iV = 16 


1.01(5) 


0.98(4) 


0.97(4) 


0.96(4) 


0.95(4) 


0.95(4) 


0.94(5) 


= 16, sme 


0.99(4) 


0.97(4) 


0.96(4) 


0.95(4) 


0.95(4) 


0.94(4) 


0.94(5) 


p (GeV) 


0.931 


1.83 


2.65 


3.37 


3.97 


4.41 


4.68 


N = 4 


- 


- 


- 


1.14(4) 


- 


- 


- 


N = A, sme 


- 


- 


- 


1.12(4) 


- 


- 


- 


N = 8 


- 


1.13(6) 


- 


1.05(5) 


- 


1.03(5) 


- 


N = 8, sme 


- 


1.03(5) 


- 


0.97(5) 


- 


0.95(5) 


- 


AT = 16 


1.04(6) 


1.00(4) 


0.99(5) 


0.99(6) 


0.98(5) 


0.98(5) 


0.99(5) 


= 16, sme 


1.01(5) 


0.97(5) 


0.96(5) 


0.95(4) 


0.94(4) 


0.94(5) 


0.94(5) 


p (GeV) 


1.29 


2.53 


3.67 


4.67 


5.49 


6.10 


6.48 


A^ = 4 








1.21(4) 








A'' = 4, sme 








1.16(4) 








N = 8 




1.17(7) 




1.05(5) 




1.03(5) 




N = 8, sme 




1.17(7) 




1.04(4) 




1.02(5) 




AT = 16 


1.04(5) 


1.00(5) 


1.00(5) 


1.00(5) 


1.00(5) 


0.99(5) 


1.00(5) 


N = 16, sme 


0.97(5) 


0.94(5) 


0.95(5) 


0.94(6) 


0.94(5) 


0.94(5) 


0.94(5) 



Table 5: Results for Z^'^{p) as a function of p — p/a (in GeV) for the three lattice sides A'' = 4, 8, 16 
and the three f3 values considered (i.e. 2.2, 2.3 and 2.4, respectively top, center and bottom part of 
the table) in the case of symmetric momenta, considering the two different gauge-fixing methods 
(without and with smearing). Error bars were obtained using the bootstrap method with 250 
samples. 





asymmetric momenta 


symmetric momenta 


/? = 2.2 


1.148(4) 


0.960(6) 


(3 = 2.2, sme 


1.056(4) 


0.952(5) 


/3 = 2.3 


1.039(3) 


0.989(3) 


(3 = 2.3, sme 


1.070(3) 


0.949(5) 


/3 = 2.4 


1.082(4) 


1.004(6) 


/? = 2.4, sme 


1.050(3) 


0.946(4) 



Table 6: Constant fits of the Z-^^ ^(p) data reported in Tables 4 and 5 for the lattice volume V = 16* 
considering momenta p>l GeV. 
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V = 4^ 


V = 8^ 


V = 16"^ 




/3 


= 2.2 


5.33(6) 


11.6(3) 


14.5(4) 






2.2, sine 


5.38(6) 


11.8(2) 


14.4(3) 






= 2.3 


6.17(7) 


18.9(4) 


31.7(7) 






2.3, sme 


6.23(7) 


19.3(4) 


32.4(7) 






= 2.4 


6.87(8) 


23.5(5) 


53(1) 






2.4, sme 


6.94(8) 


23.8(5) 


53(1) 



Table 7: Results for the gluon propagator at zero momentum D{0) (in lattice units) for the three 
lattice sides (TV = 4, 8, 16) and the three (3 values (i.e. 2.2, 2.3 and 2.4) considering the two different 
gauge-fixing methods (without and with smearing). Error bars were obtained using the bootstrap 
method with 250 samples. 





V= I' 


V = 8i 


\' = 16'^ 


13= 2.2 
P = 2.2, sme 


0.353(1) 
0.352(1) 


0.3174(5) 
0.3175(5) 


0.3098(3) 
0.3097(3) 


/3 = 2.3 
P = 2.3, sme 


0.340(1) 
0.340(1) 


0.3075(4) 
0.3075(4) 


0.2998(2) 
0.2997(2) 


f3 = 2.4 
P = 2.4, sme 


0.334(1) 
0.333(1) 


0.3013(4) 
0.3013(4) 


0.2946(3) 
0.2945(3) 



Table 8: Results for the ghost propagator G{p) (in lattice units), evaluated at the asymmetric 
momentum with components pi = P2 = Ps = and P4 = N/2, for the three lattice sides {N = 
4,8,16) and the three (3 values (i.e. 2.2,2.3 and 2.4) considering the two different gauge-fixing 
methods (without and with smearing). Error bars were obtained using the bootstrap method with 
250 samples. 



P (CeX) 


n..H66 


0.718 


1.04 


1.33 


1.56 


1.73 


1.84 


TV = 16 
N = 16, sme 


0.28(2) 
0.26(2) 


0.019(1) 
0.0177(9) 


0.0041(2) 
0.0038(2) 


0.00143(7) 
0.00131(7) 


0.00062(3) 
0.00057(3) 


0.00029(2) 
0.00027(1) 


0.000122(6) 
0.000111(6) 


p (GeV) 


0.466 


0.913 


1.33 


1.69 


1.98 


2.20 


2.34 


iV = 16 
N = 16, sme 


0.46(3) 
0.48(4) 


0.031(2) 
0.032(2) 


0.0069(4) 
0.0073(4) 


0.0025(1) 
0.0026(1) 


0.00113(6) 
0.00118(7) 


0.00054(3) 
0.00057(3) 


0.00023(1) 
0.00024(1) 


p (GeV) 


0.644 


1.26 


1.83 


2.34 


2.75 


3.05 


3.24 


N = ie 

N = 16, sme 


0.74(7) 
0.70(6) 


0.046(3) 
0.045(3) 


0.0108(6) 
0.0106(6) 


0.0040(2) 
0.0039(2) 


0.0018(1) 
0.0018(1) 


0.00089(5) 
0.00087(5) 


0.00037(2) 
0.00037(2) 



Table 9: Results for the quantity J2a,b,c (^m(^) (M'^Y" (p) ) as a function of p = p/a (in 

GeV) for the lattice side N = 16 and the three /3 values considered (i.e. 2.2, 2.3 and 2.4, respectively 
top, center and bottom part of the tabic) in the case of asymmetric momenta, considering the two 
different gauge-fixing methods (without and with smearing). Error bars were obtained using the 
bootstrap method with 250 samples. 



- 18 - 



